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Abstract. In this article, results have been presented for the two-time correlation functions for a free and a harmon- 
ically confined Brownian particle in a simple shear flow. For a free Brownian particle, the motion along the direction 
of shear exhibit two distinct dynamics, with the mean-square-displacement being diffusive at short times while at late 
times scales as ? 3 . In contrast the cross-correlation (x(t)y(t) ) scales quadratically for all times. In the case of a harmon- 
ically trapped Brownian particle, the mean-square-displacement exhibits a plateau determined by the strength of the 
confinement and the shear. Further, the analysis is extended to a chain of Brownian particles interacting via a harmonic 
and a bending potential. Finally, the persistence probability is constructed from the two-time correlation functions. 
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1 Introduction 

The phenomenon of persistence has been well studied over the 
past decade, both theoretically and experimentally |U]|2]. Per- 
sistence is the probability that a stochastic variable x(t) retains 
a particular property up to the observation time interval t - the 
property being either the sign of the variable or the crossing of 
the origin J2][3j|4l. The later definition of p(t) is synonymous 
with the survival probability of the stochastic process, and the 
problem can also be approached using the backward Fokker- 
Planck equation. For a wide class of non-equilibrium systems, 
the asymptotic decay of p(t) exhibits a power-law decay with 
a non-trivial exponent 9. This algebraic decay, and the expo- 
nent, has been investigated for a wide class of non-interacting 
as well as interacting systems including the overdamped Brow- 
nian walker in an infinite [5 1] and finite medium [6], diffusion 
equation with random initial conditions 010 > advection of a 
passive scalar [8], fluctuating interfaces ll WT0l[TTll . critical dy- 
namics iflZl . granular media lfT3l[T4l . disordered environments 
lll5lll6lfT7l and polymer dynamics [4J. The probability p{t) 
acts as a dynamic probe, indicating how the system retains the 
memory of its initial configuration as it evolves in time, and is 
also an indirect test for the two-time correlation functions for a 
non-stationary process. However, the estimation of the persis- 
tence probability is notoriously difficult and an exact analyti- 
cal prediction for p(t) at all times can be made only when the 
stochastic process is Gaussian and Markovian, such as, for an 
overdamped Brownian particle. 



When the stochastic dynamics is a Gaussian Markovian 
process, the non-stationary process x(t) can be mapped to a sta- 
tionary Ornstein-Uhlenbeck process X(T), for which, the sta- 
tionary correlator C(T) = (X (T)X (0) ) is exponentially decay- 
ing at all times. Asymptotically, the survival probability of the 
stationary Ornstein-Uhlenbeck process is proportional to C(T) 
and p(t) can be obtained by the inverse time transformation 
applied to X lfTl [T8l . For a wide class of systems, the station- 
ary correlator C(T) is often non-exponential and this straight- 
forward method can no longer be applied. Using the classifica- 
tion of Slepian ifTSl . a correlator is of class a, if C(T) ~ 1 — 
ff(T a ) in the limit of T — >• 0. The behavior of C(T) near zero 
characterizes the density of zero crossings p of the stochastic 
process X [318 19|. When a = 2, the number of zero cross- 
ings is finite and the Independent Interval Approximation can 
be used to predict the exponent 9 [3|. If however, a < 2, the 
stochastic process has infinte number of zero crossings and a 
perturbative expansion gives a fairly good estimate for 9 ||9l. 

The case of a free Brownian particle, in the overdamped 
limit, is particularly simple and illustrates how the persistence 
probability can be determined for a Gaussian Markov process 
fl8l . We define the persistence probability as the probability 
that the particle has not crossed the origin during the observa- 
tion time interval t . Even though the system under considera- 
tion is very simple, its application is abundant in quantitative 
science [20|. The two-time correlation function for the non- 
stationary process x{t) = Jq T7(f')d?', where 77 is a Gaussian 
stochastic noise, can be transformed into a stationary Ornstein- 
Uhlenbeck process using the successive transformations X = 
x/ \J (x 2 (t) ) and T = Inf. The stationary correlator of X is 
then given by e~ T l 2 , and following Ref. lfT8l . the persistence 
probability in the transformed variable decays exponentially - 
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P{T) ~ e- T < 2 . Using the reverse transformation, the decay of 
p(t) in real time becomes algebraic with an exponent = 1/2. 

In this article, we investigate the persistence probability of 
a free and confined Brownian particle in a shear flow and that 
of an interacting chain of Brownian particles. The work is mo- 
tivated on one hand by the extensive theoretical [21 ,22,23,24, 
25 1 and experimental [26 1 study of Brownian motion in a linear 
shear flow and their importance in micro fluidic applications, 
and on the other hand, due to the rich dynamics exhibited by 
such systems as a result of an interplay of thermal fluctuations 
and the imposed velocity gradient [27 28 29]. The investiga- 
tion of the survival probability for a free Brownian particle in 
a deterministic flow field v(x) has been previously investigated 
[30 ,31 32 1, and is known that for all odd functions of v(x) the 
survival probability decays as f~'/ 4 . However, all of these ap- 
proaches use the Fokker-Planck equation to determine the sur- 
vival probability. In the present work, we emphasize on the 
two-time correlation function and present explicit results for 
them for the case of a free and harmonically confined Brown- 
ian particle in a transverse flow field. In the case of a chain of 
Brownian particles harmonically bounded to its nearest neigh- 
bors, we impose an additional bending potential to mimic the 
case of a semi-flexible polymer chain. 

The rest of the article is organized as follows: the model 
system is introduced in section [2] The relevant results for the 
two-time correlation functions and the persistence probability 
for a free Brownian particle is presented in section [2] and for a 
confined Brownian particle, in section [3] We extend the analy- 
sis presented in the previous sections to a chain of interacting 
Brownian particles and present the results in section|4] 



2 Brownian particle in a shear flow 

The stochastic dynamics of a Brownian particle can be looked 
at from different levels of coarse-graining. Typically, if the mea- 
surement time intervals of the position and momenta of a col- 
loid are well separated, the Markovian Langevin equation in 
the overdamped limit faithfully reproduces the dynamics. In 
reality, however, the dynamics of a Brownian particle exhibit a 
far more rich physics. As opposed to the Markovian Langevin 
equation, characterized by a single relaxation time, the actual 
dynamics of Brownian particle is characterized by a set of time- 
scales determined by sound propagation, vorticity diffusion and 
the inertia of the colloid. As a consequence, the motion of the 
Brownian particle no longer remains Markovian and the gener- 
alized Markovian equation with a memory dependent friction 
and a correlated noise is used to describe the dynamics. With 
this increasing level of complexity, an analytically tractable re- 
sult becomes difficult and one has to resort to a numerical inte- 
gration of the equations of motion. In the following, however, 
we consider the overdamped Langevin equation to describe the 
dynamics, neglecting the effects of solvent and the inertia of 
the Brownian particle. 

We consider the motion of a Brownian particle, with a unit 
mass, in an unbounded solvent moving in a two-dimensional 
planar geometry. Following [JQ, we consider a stationary dis- 
tribution of velocity, 

u=(0,fljt). (1) 



The force on the Brownian particle, due to the imposed flow is 
given by F = — £(v — u), where v is the instantaneous velocity 
of the particle and £ is the Stokes's friction on the colloid. The 
Langevin equation for the position of a colloid r = (x.y), in the 
overdamped limit, takes the form 

— = rk(t) and ^ =ax+r]y(t), (2) 

with 77 = (t] x , rjy) as a Gaussian white noise with correlations 

(77(f)) =0 and (i]{t)®T){t'))=2Dl8{t-t'), 

(3) 

where I is the identity matrix and ® denote the outer product of 
a vector quantity. The strength of the noise correlations is given 
by the diffusion constant D = IcqT /£. The two-time correlation 
functions from equation |2) becomes, 

(x(h)x(t 2 )) = f^t[ J\t' 2 ( nx {t[) nx {t' 2 )) (4) 

and 

(y(h)y(ti)) = / o V J\t' 2 (riyit'Mti)) 

+ a 2 f l dt[ f' 2 dt' 2 (x(t[)x{t' 2 )) . (5) 
Jo Jo 

Assuming t\ > t 2 , equation |4} yields (x(t\)x(t 2 )) = 2Dt 2 and 

(y(h)y(t 2 )) =2Dt 2 + a 2 [ l 6t[ [ 2 dt 2 mm{t[,t' 2 ). (6) 
Jo Jo 

The integral on the right-hand-side of equation |6} evaluates to 
(y(ti)y(t 2 )) = 2Df 2 + a 2 Z)(vf-|) . (7) 

Similarly, the cross-correlation functions (x{t\)y(t 2 ) ) and (y{t\)x(t 2 ) ) 
can be constructed as, 

(x(h)y(t 2 ) )=a£ dt[ £ dt> 2 < J] x (t[)x(t 2 ) ) 

+ f'drj [' 2 dt> 2 (riMHVi)) , (8) 
Jo Jo 

and 

{y(h)x(t 2 ))=aj 1 dt[j Z dt' 2 (x(t[)r] x (t' 2 )) 

+ J\[J\t 2 {n x {t[)T)y{t' 2 )). (9) 

The second terms in both of the above equation are zero due 
to the choice of the noise correlations in equation d3j. Further, 
using the solution x(t) = /odf'rj^?'), the correlation function 
in equation <(8j becomes, 

(x(h)y(t 2 )) =a £dt' 2 pdt[ J\t' 2 ' (r tx {t[)7]{t , 2 , ))=Datl 

(10) 
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and the correlation function in equation (|9]l evaluates to, 

(y(t 1 )x(t 2 ))=a f l dt[ f' 2 dt' 2 /V/ (Tkffiikti)) 
Jo Jo Jo 



Using the time transformation e T = t, equation (jT3J and equa- 
tion ( |T6] > are transformed to a Gaussian stationary process with 
correlations, 



2Da j dt[ min(f2,f' 1 ) = 2Da ( t\t 2 - J- 



C(T) 



e -T/2 

l e -T/2_l_ e -3T/2 



for t <C a 1 
for t > cT l . 



(17) 



(11) 



Using the two-time correlation functions evaluated above, 
we can construct the corresponding mean-square displacements 
by simply substituting t 2 = t \ = t. Along the x-direction the mo- 
tion remains purely diffusive, while the mean-square-displacemen 
of the particle along the y-direction becomes, 



(y 2 (t))=2Dt + -a 2 Dt 3 . 



(12) 



For short times, t << a -1 , the motion of the Brownian particle 
along the y-direction is purely diffusive, while in the asymp- 
totic regime of t >> a -1 the mean-square-displacement scales 
as f 3 - similar to that of a randomly accelerated particle. 

The cross-correlation function (x(t)y(t) ), however, has two 
contributions and adding up equation (jTOJ and equation ( fTT| 
gives, 

(x(t)y(t)) = Dat 2 . (13) 



Our aim is to extract the persistence probability of the Brow- 
nian particle using the correlation function defined in equa- 
tion Q. To this end, we define it as the probability that the 
y-component of the particle position has not crossed zero. The 
fundamental idea is to map the non-stationary process y(t ) to 
a stationary Ornstein-Uhlenbeck (O-U) process Y. The station- 
ary correlator C(T) for Y decays exponentially for all times and 
the persistence probability can be shown to decay as P(T) ~ 
| sin -1 [C(J)] Q8 ,33 1. The mapping to the stationary O-U pro- 
cess is done using two transformations, first normalizing the 
stochastic variable by the root-mean-square displacement and 
then a suitable transformation in time. For the corr elation f unc- 
tion in equation (FT), the transformation Y = y(t )/ \J (y 2 (t) ) yields 



Since the stationary correlator for t < a 1 is exponentially de- 
caying, the persistence probability in the transformed variable 
T is given by P(T) ~ e~ T l 2 and in real time p(t) ~ f~'/ 2 . In the 
asymptotic regime, the correlation function corresponds to that 
of a randomly accelerated particle and the persistence proba- 
bility is known from the works of Sinai and Burkhardt B34II141 
to decay as pit) ~ t _1 / 4 . 

To validate the results presented above, numerical simu- 
lations were done by integrating equation |2} using the Eu- 
ler scheme with an integration time step of df = 0.001. From 
this numerical integration, the position of the particle in the 
subsequent time steps were measured and the mean-square- 
displacement along and perpendicular to the direction of the 
applied shear, as well as the cross-correlation functions were 
computed. The results of the numerical simulation is presented 
in the left panel of fig.[T]and compared to the analytical results 
of equation (jT2]» and equation ( |13) , The measured data were 
averaged over 10 5 independent configurations. In all of these 
simulation runs, the particle was started in the neighborhood of 
zero, so that the sign of y(Q) is defined and subsequently the 
sign of y(t) was followed. The fraction of particles for which 
the sign of y(t) did not change quantified the persistence prob- 
ability. The measured persistence probability for the Brownian 
particle in the direction of the applied shear is shown in the 
right panel of fig.[T] The measured probability exhibits two dis- 
tinct regimes of decay, for t << a -1 a decay of f -1 / 2 and for 
t » cT x a decay of t 



2Dt 2 +a 2 D 



(Y(h)Y(t2)) 



ht 2 - 



[2DU 



% 2 Dt~ 



,1/2 



(2Dh 



(14) 

A suitable time transformation which converts equation ( fT4} 
to a stationary process for all times is non-trivial. However, 
we can extract persistence probability by considering the lim- 
iting behaviors of equation ( |14) , As noted above, for t <C a -1 , 
the motion of the colloid is purely diffusive and neglecting the 
<??(f 3 ) terms gives, 



(Y(h)Y(t 2 )) 



(15) 



In the opposite limit of t 3> a -1 , the linear term can be ne- 
glected and the correlation function for the normalized variable 
becomes, 



(Y( tl )Y(t 2 )) = 



3 ft 2 



2 Wi 



1/2 



1 (t 2 W 

2 U 



(16) 



-1/4 



3 Harmonically confined Brownian particle in 
a shear flow 

In this section we study the dynamics of a tracer in a sim- 
ple shear flow which is also confined by a harmonic potential. 
The harmonic confinement occurs naturally in the experiments 
when a tracking of a tracer is done using an optical tweezers. 
The equations of motion for the position of a colloid which is 
subjected to a harmonic confinement U (r) = jKr 2 and the sta- 
tionary velocity profile of equation ([T} are, 

— = -Kx + r\ x (t) and -j^ = -Ky + ax + T\ y {t), (18) 

together with the noise correlations of equation Q. The time 
evolution of the coordinates is given by, 



x(t) = / df V 
Jo 



t} Vx{t') and 
y (t) = ^6t'e-< t -^[x(t')+T 1 y(t')]. 



(19) 



Assuming that t\ > t 2 , the two-time correlation function for the 
variable x(t) is straightforward and gives, 

D r 



(x{h)x(t 2 )) 



e *(ti-t 2 ) _ e -K(t l +t 2 ) 



(20) 
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Fig. 1. (a) Mean-square displacement of a Brownian particle in a simple shear along the y direction for a = 1.00 (•) and 5.00 (■). The 
corresponding mean-square displacement along the ^-direction is shown for a = 1.00 (O) only. The solid lines are the plots of equation i ] 1 2^ 
for the respective values of a and the dot-dashed line is the plot of (x 2 (?) ) = 2Dt. Inset: Plot of (x(t)y(t) ) for a = 1.00 ( A) and 5.00 (V). The 
dashed line are the plots of equation < | 1 3| > . (b) Persistence probability of a Brownian particle in a simple shear for (•), 5.00 (■) and 10.0 (A). 
The dashed line is plot of ?~ 1//2 corresponding to that of a free Brownian particle, while the dot dashed line is a plot of f -1 / 4 . 



while fory(f) we have 
a 2 D 



(y(ti)y(hj 



e «H+tz) r' i dt < f 2 & e K(,' l+ ,' 2 

K Jo Jo 

(21) 



p -K{t[+t> 2 ) 



The integral over t [ and t' 2 yields, 



(y(fi)y(fe)) = 



a 2 D 
2>? 



,-<h-t 2 ) _ -K(f I +r 2 ) 



K((h+t 2 )e-"^-( tl -t 2 ) 

-K{t\-t 2 ) _ e -K{fl+t%) 



- 2k 1 tyt 2 e-< h+ ^ 



D 

K 



(22) 



The mean-square displacement for y(t) immediately follows 
from equation p2| ) via a substitution t\=t 2 = t, 



(y 2 (t)) = ^ [{\-e- lKt )+2Kte- la -iK^e- 1 *^ 



--\l-e- 2Kt ] 



(23) 



Further, the cross-correlation functions (x(ti)y(t 2 )) and 
(y(h)x(t 2 )) is given by, 

<*(*iMfe)>= / o ' 2 d^ e-«^ {x{h)x{t> 2 )) 

= a® \ e -«h-t 2 ) _ e -<t 1+t2 ) _ 2Kt2e - K ( tl +t 2 ) 

2k 2 L 

(y(h)x(t 2 )) = £dt[ e-«i-V (x(t[)x(t 2 )) 



(24) 



aD r f i 

K Jo 



(25) 




Fig. 2. Mean-square displacement of a Brownian particle in a simple 
shear along the y direction for a = 1.00 (•), 2.50 (■) and 5.00 (A). 
The corresponding mean-square displacement along the x-direction is 
shown for a = 1.00 (O) only. The solid lines are the plots of equa- 
tion {23} for the respective values of a and the dot-dashed line is the 
plot of(x 2 (t) ) = 2Dt. Inset: Plot of (x(t)y(t) ) for a = 2.50 ( A) and 
5.00 (~). The dashed line are the plots of equation J27l. 



The integral over t[ in equation ( |25] > gives 
aD 



(y(h)x{t 2 )) 



2k 2 



-<h-t 2 )_ e -K{ h+ t 2 ) +lK ( h - t2 ) e -<h-t 2 ) 



-2Kt\e 



-K{h+t 2 ) 



(26) 



The cross-correlation function (x(t)y(t) ) is obtained by adding 
equation p4| and equation (p6]l and the substitution t\=t 2 =t 
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Fig. 3. Artistic illustration of a bead-spring model of a polymer. In 
addition to the harmonic confinement between the beads, an angular 
confinement is also implemented. 



which gives, 



?- 2w -2Kte- 2Kt ] 



(27) 



A Taylor expansion of equation ( |23| l for t < fc _1 , shows 
that the dynamics at short times scales as 2Dt + (2/3)a 2 Dt 3 , 
while in the asymptotic regime, the mean-square displacement 
saturates to a value D/K + a 2 D/K 3 . To determine the persis- 
tence probability, we follow the steps outlined in the previous 
section. However, an expression for the persistence probabil- 
ity valid at all times is once again not feasible. Even in the 
asymptotic regime, a suitable time transformation does not ex- 
ist which converts the process y(t) to a Gaussian stationary pro- 
cess. Progress can only be made in the time domain which is 
smaller than K^ 1 , where the two-time correlation function is 
identical to equation (IT). Accordingly, the persistence prob- 
ability shows an initial decay of f~ 1//2 , followed by a decay 
of f -1 / 4 . The mean-square-displacement of the harmonically 
confined Brownian particle in a simple shear flow is shown in 
fig. [2] Once again, the integration of equation ( fT8| i was carried 
out using the Euler scheme with an integration time step of 
d? = 0.001. The measured mean-square-displacement is com- 
pared to the analytical prediction of equation ( |2"3j ) in the main 
figure, while the inset depicts the cross-correlation function 
(x(t)y(t) ) for two different values of a, defined in equation ([l}. 



4 Interacting chain of Brownian particles 

In this section, we investigate the persistence probability of an 
interacting chain of N Brownian particles, in a two-dimensional 
planar solvent subjected to the stationary flow of equation ([TJ. 
In addition to the harmonic confinement U(r) = ^K(r — Rq) 2 
of a particle to its nearest neighbors, we impose a bending po- 
tential 



U 



bend 



N-l 



(0)=£^(l+u,u ; --i) 2 , 



(28) 



where u, = r 1+ i — r,-/|r,-+i — r,-|, This discrete model of a chain 
of Brownian particle is often used in simulations to model semi- 
flexible polymer chains [35 .36 1. In the absence of the bending 
potential, the model corresponds to that of a Rouse chain and 



the persistence probability of such a chain in a simple shear 
flow has been studied by Bhattacharya et. ai, 0). The resulting 
expression for the bending force takes the form, 



f 



bend 



Kg £(l+ft^l)3r¥j 
7=2 ° Ti 



-1 



(29) 



Using the explicit form of u,, and assuming that the confine- 
ment strength of the harmonic force and the bending rigidity is 
sufficiently large so that we are in weakly bending limit, such 
that |u(z + 1) — u(/)| <C 1 [37 1, the force on particle i becomes, 



f 



bend = _ 2 g ( _ r ._ 2+4r ._ 1 _6r,+4r i+1 -r /+2 ) 
"o 



(30) 



with Ro as the inter-particle separation. Note that in this ap- 
proximation, the force due to the harmonic confinement can be 
neglected. The equation of a Brownian particle then becomes, 

-2 ^j(-*/_2 + Ax i- 1 ~ + 4x m ~ x i+2) + ^xiht) 
Ri 



o 



dt 

= -2-r(-yi-2 + 4yi-i-6yi+4yi + i -yi+i) +axi(t) + J? y (/,f) 

(31) 



dt 



The above equations are only valid for N — 4 particles, and the 
4 particles at the boundary will have different equations of mo- 
tion. However, for an infinitely long chain, the monomer dy- 
namics at late times is not sensitive to the boundary conditions 
and accordingly a continuum version of the equations can be 
formulated by replacing i with a continuous variable s, 

d 4 x(s,t) 



dx(sJ) 

— V* 

dt ds 4 



(32) 



d 4 y(s,t) 



-ax(t) + rj y (s,t) 



dy{s,t) = _ w 
dt ds 4 
Fourier transforming the above equations leads to, 

dx 

dy 
dt 

We have dropped the stochastic noise term in equation ( j34] >, 
since at late-times the second term ax(t) becomes the dominant 
term compared to the stochastic force. The formal solutions of 
equation (|33]l and equation 04} are given by, 



-Kk 4 x(k,t) + i) x (k,t) 
-K k 4 y{k,t) +ax(k,t) 



(33) 
(34) 



x(k,t) 



e -Kk* {t-t>) 



y(k,t)= f dt' e- Kk4 ('-'') x(k,t'). 
Jo 



(35) 
(36) 



Using the noise correlation for (f\ x (k 1 t)f] x (k l ,t)) = 2D8(k- 
k' ) 8 (t — t ' ) , the two-time correlation function for x reads, 



(x(k,h)x(k' ,t 2 )) 



D 

W) 

D 

W) 



8{k + k')e^ k){ - h+tl ' ) [e^WH'iA) _ \ 



8{k + k') 



-D{k){t x -t 2 ) . 



-D{k){h+t 2 ) 

(37) 
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where D(k) — Kk 4 is the diffusion coefficient of a mode k. The 
second line of equation ( |37| ) assumes fi >t 2 . Similarly, the two- 
time correlation function for the y(s,t) can be constructed from 
equation ( |37j ), 



(y(s,tx)yM) = 



{ki+k 2 )s 



dt\ / dt' 7 



dki f dk 2 

2k J 2k jo Jo 

e - 5 (* 1 )(' 1 -''l) e - 5 fe)fe-4)(y( fel) ^)y( fe) 4)) (38) 

Substituting equation p7) > in the equation above and subse- 
quently integrating over the delta function and renaming the 
dummy variable k\ we arrive at, 



(y(s,h)y(s,t 2 )) =af 
D 



dk f'i 



2k Jo 



dt\ 



■ ^ e -D{k){ h +t 2 ) 



D(k) 



a z D f'i 



dt\ / dt' 7 



-1 

dk l- e -^ 4 ('i+'2) 



2k 



k 4 



dk 1 - e -' dr4 ('i+'2- 2min (''i,f2)) 



2k 



k 4 



(39) 



The integral over k yields, 

,2 



(y(s,h)y(s,t 2 )} 



K Jo JO 



(fl+?2) i 



-(ti+t 2 -2mm(t[,t' 2 )) i 



The integral over the t [ and t' 2 finally yields, 
2a 2 Dr(i) 



(y(s,h)y(s,t 2 )} 



3k 



tMh+t 2 )i - —(h 

2 , ,n 
+ TT (fi-f 2 ) T 



(40) 



-h )t 



(41) 



Note that the correlation function is independent of s, and this 
translational invariance is a consequence of taking the limit of 
an infinitely long chain. We shall, henceforth, drop the con- 
tinuous variable from the notations. At this point, we justify 
dropping the stochastic noise term in equation ( (34] i. The inclu- 
sion of this term leads to an additional term ^(V 3 / 4 ), which can 
be neglected compared to the terms in correlation functions in 
equation (|4l}, which are of G{t n l 4 ). 

The transformation to a stationary process is done using the 
identical procedure as outlined in the previous sections. Defin- 
ing the normalized variable Y — y(t)/y/ (y 2 (t)) and the trans- 
formation e T — t yields a stationary correlation function, 



not applied here. A Taylor expansion of C(T) in the neighbor- 
hood of zero gives, 



C(T)=(Y(0)Y(T)) = 



11 






* 8 






cosh | 


}): 




cosh [ — | 


y 






~3 


w. 



T 
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(42) 



Since the correlation function is not exponentially decaying, 
the simple route to determine the persistence probability can 



C(T) = 1 



55 2 
96 



■0{T n l 4 ), 



(43) 



demonstrating that Y is a smooth process with a finite num- 
ber of zero crossings. The density of zero crossings is given 

|| and the mean interval 



by mm, p = v^cW)/* = \ 

size (T) is given by (T) = 1/p. To determine the exponent 
9 we use the Independent Interval approximation [3|, where 
the intervals between successive zero crossings are assumed to 
be independent. The distribution P(T) of the intervals between 
the successive zeros can be related to the stationary correlator 
C(T) in the Laplace domain as JH, 



P(s) 



1 



(44) 



where the tildes refer to the Laplace transforms of the corre- 
sponding quantities. Finally, the estimation of the exponent 9 
from the interval size distribution translates to the determina- 
tion of the simple pole of P(s), that is, the root of the denomi- 
nator in equation (j44j» |4][3], 



F(s) = 1 + 



'J)s 
2 



s I dTe 
lo 



-sT 



cosh 



sinh 



cosh 



(45) 



The numerical estimation of the root of F(s) yields a value 
of 9 = 0.29695. To verify this result, we performed numerical 
integration of equation (J32j by discretizing the equation on a 
lattice with 1024 lattice points, 

Xi (t m+ l)=Xi (t m ) — AtK [-Xi- 2 (t m ) + 4x;_ i (t m ) - 6x t (t m ) 



+4x;+i (f m ) -x i+ z(t m )] + V2DAtJ^(t m ) 

yi (tm+ 1 ) = yi (tm :) - ^ t K [-)?/_ 2 (t m ) + 4y i (t m ) - 6yi (t m ) 

)}+axi(tm), (46) 

where JYiifrn) is a Gaussian random number with zero mean 
and unit variance and At = 0.01 is the integration time step. 
The free boundary conditions were implemented by choosing 
r(0) = r(l) = r(2) and r(N) =r(N-l) = r(N-2). In the 
course of the simulation, the sign of the variable y,(f) were 
monitored and the persistence probability was defined as the 
fraction of the lattice points for which the sign of yi(t m ) was 
same as that of y, (0). Further, the measured persistence proba- 
bility was averaged over 500 independent configurations. The 
asymptotic of the probability decays as a power law with an 
exponent close to the prediction of IIA (see fig.|4|. 

5 Conclusion 

In conclusion, we have presented results for the two-time cor- 
relation functions of a single free and confined Brownian par- 
ticle in a simple shear flow. The confinement of the Brownian 
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Fig. 4. Persistence probability of a semi-flexible polymer in a shear 
flow for values of a = 0.1 (•) and 0.50 (□). The solid black line is 
a plot of with = 0.29695, the prediction of II A. A power law 
fit to the asymptotic of the data yields 6 = 0.306236 for a = 0.1 and 
= 0.297205 for a = 0.5. 



particle is modeled as a harmonic trap, often encountered in 
trapping and tracking experiments. The persistence probabil- 
ity, defined as the probability that the sign of the stochastic ob- 
servable has not changed sign up to time t is constructed from 
the correlation function. The probability shows two distinct al- 
gebraic decays. For short times t « a -1 , when the particle 
does not feel the effect of the shear, the motion is found to be 
purely diffusive while at late times the motion is super ballistic 
with the mean-square-displacement along the direction of shear 
scaling as f 3 . We have also extended the analysis to a chain of 
Brownian particles interacting via a harmonic potential and a 
bending potential. The asymptotic of the persistence probabil- 
ity is found to decay as a power law with as exponent in close 
agreement with that of the predictions from Independent Inter- 
val Approximations. 
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